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Abstract 

We study the stability of the rephca-symmetric solution of a two-sublattice infinite-range spin- 
glass model, which can describe the transition from antiferromagnetic to spin glass state. The 
eigenvalues associated with replica-symmetric perturbations are in general complex. The natural 
generalization of the usual stability condition is to require the real part of these eigenvalues to be 
positive. The necessary and sufficient conditions for all the roots of the secular equation to have 
positive real parts is given by the Hurwitz criterion. The generalized stability condition allows a 
consistent analysis of the phase diagram within the replica-symmetric approximation. 
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I. INTRODUCTION 



The infinite-range Sherrington-Kirkpatrick fSK) model jl| for a spin glass has attracted 



These investigations have revealed 



considerable attention over the past decades 2, _, 
highly non-trivial properties such as the instability of replica-symmetric fRS}solution and 
the replica-symmetry-breaking scheme to produce a stable solution j^I?,!^ ^- Most studies 
have concentrated on situations where the exchange distributions are either symmetric or 
with an additional ferromagnetic interaction. More recently a two-sublattice version of 



the SK model was introduced 



to allow for antiferromagnetic interactions 



between different sublattices. Such extension is quite natural in view of the existence of 



many experimental systems such as Fe2:Mgi_a,Cl2 H Q, Q and Fe,.Mni_,Ti03! 
which exhibit transition from and Ising antiferromagnetic into an Ising spin glass state for 
certain range of x values. In contrast to the standard SK model, in the two-sublattice SK 
model with antiferromagnetic intersublattice interactions, the ordered (antiferromagnetic) 



17 



phase extends to finite fields and the de Almeida-Thouless instability line 



has distinct 



branches in the paramagnetic and antiferromagnetic phases, which do not meet at a first- 



order transition 



|l2,ll3|. Experimental determination of the field-temperature phase 

are in 



diagram in Fea:Mni_a;Ti03, as well as the de Almeida-Thouless instability line 
qualitative agreement with mean- field results [iS^ . 

In the previous studies of this model the stability of the RS solutio n ag ainst transversal 



fluctuations, i.e., outside the RS space, has already been investigated (H 
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and 



the stability against longitudinal fluctuations, i.e., inside the RS space, was also briefly con- 



sidered 



12| . The stability of the RS solution against transversal fluctuations is important to 



establish whether replica symmetry breaking is necessary. The stability against longitudinal 
fluctuations, however, is also necessary to ensure the validity of RS solution. For certain 
parameter values of the two-sublattice SK model there may be up to three RS solutions, 
all of them stable against transversal fluctuations. In such a situation the analysis of the 
stability against longitudinal fluctuations is important for a consistent study of the phase 
diagram by eliminating unstable solutions. 

In this work we remedy the lack of such investigation by a detailed numerical and ana- 
lytical study of the eigenvalues associated with longitudinal fluctuations. Surprisingly, these 
eigenvalues are in general complex. It is natural to assume that stability condition should 
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require the real part of these eigenvalues to be positive. The necessary and sufficient condi- 
tion for all the roots of the secular equation to have positive real part is given by the Hurwitz 
criterion. We show that this generalized stabihty condition allows a consistent study of the 
phase diagram within the RS approximation. 



II. THE MODEL 



We consider a system of 2N Ising spins Si = ±1 located at the sites of two identical 
sublattices A and B. The interactions are described by the Hamiltonian 

'hi = JijSiSj J^jSiSj J^jSiSj H Sj, (1) 

ieA,j€B («i)eA fo')GB i 

where the ffist sum is over all distinct pairs of spins belonging to different sublattices, the 
second and third ones refer to all distinct pairs of spins belonging to the same sublattices, and 
the last sum is over all spins in the two sublattices. Jij is the exchange interaction between 
spins in different sublattices, J-j is the exchange interaction between spins in the same 
sublattice, and H is the applied magnetic field. The exchange interactions are independent, 
quenched, Gaussian random variables with mean values 

Jo J' 

{Jij) J = {J'ij)j = (2) 

and variances 

The mean intrasublattice interactions will always assumed to be ferromagnetic (Jq > 0), 
whereas the mean intersublattice interactions may be ferromagnetic (Jq > 0) or antiferro- 
magnetic (Jq < 0). 

The standard approach to compute the quenched average is to introduce n non-interacting 
replicas a = l,2,...,nof the system, calculate the annealed averages and then take the limit 



Urn 

In this replica method the free energy per spin / is given by 

/ = ^ lim -0, = - lim — In (Z") , , (4) 

where (3 = l/k^T and is the partition function of n replicas of the system. Performing 
the average of over the random couplings we find 



/3Jo m'Xm'^ 



(a/3) 



2 7/2 



lP.T 



(5) 



where (a/5) denotes distinct pairs of replicas and we have introduced the sublattice magne- 
tization and sublattice overlap function of the replicas, 



m 



X 



1 



(6) 



The trace over the spin variables in ^ can be performed by taking into account the con- 
straints (jni) by means of the identities 



1 



drrix 



-too 2m 



exp 



and 



dq 



NdX 



X 



X 



2m 



exp 



-NXf qf 



N 



Ena of: 



iex 



{X = A,B), (7) 



iX = A,B). (8) 



We then obtain 

/c 



dm^. 



^ dm% ^ n / dqf / 

-ioo Zm J —oo J —too Z7C1 , r,\ J ~oc J —i' 



Nd\% 



fioo NdX 



A 



2m 



dql 



.p n^^ NdX 



a/3 
B 



2m 



where 



2 ^^""1^- 2 



exp 



n 
N 



f3J' 



-N<P{m% m% qf, qf; \% X% Xf, A^^)! , (9) 



a a 



(al3) 



(al3) 



+ E i^A^nl + Km'i) + E {^A^A + ^i^qf) - In Tr exp Ha - In Tr exp Um 



(a/3) 

with Ti^ and 1-Lb denoting the "effective sublattice Hamiltonians" 

nx = Y.^xS'" +Y. ^fS^'Sf^ {X = A,B). 

a (a/3) 



(11) 



In the limit of large the integrations over the A variables in (p)) can be performed by the 
saddle-point method. The saddle point is given by 
Tr 5" exp Hx 



m^x 



Tr expK 



X 



„fl Tr5"5^exp7^x 

qx = = • 

Tr exp7-^x 



{X = A,B). 



(12) 
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These equations determine A variables in terms of m and q variables. The remaining inte- 
grations over the m and q variables in Q can be performed by the Laplace method in the 
limit of large N. The stationary-point equations are given by 

= + J^m^ + /? Jom^, xf = (3^ J' + (3' q!^^ , {X = A,B), (13) 

where X is the sublattice complementary to X, i.e., ii X = A then X = B, and vice versa. 
Substituting these results in the expression of given by Eq. (fTUI) we find 



+ 



2 



" (a/3) 



In Tr exp Ha — In Tr exp Ti^, 



(14) 



where we have discarded terms that vanish in the limit of large N. Analogously, the effective 
sublattice Hamiltonians (fTT| become 

nx=(3j2{H + JornS, + Jom^) + E {j' '^x + S^S^ (X = A,B). (15) 

« (a/3) 

To evaluate the general expressions obtained thus far it is necessary to impose some 
structure on m and q variables. The simplest assumption corresponds to the RS solutionj^Ol 
obtained by assuming order parameters independent of replica indices, 

m^=mx, gf = gx, {X = A,B). (16) 

Proceeding in the usual way , one finds that the saddle-point equations (fT^ and 

stationary-point equations (fT^ give the equations of state 



mx = (tanhi^x) 



qx = (tanh Hx 



{X = A,B), 



(17) 



where 



Hx = f3(H+ J^nix + Jomj^+ Jj'^qx + J^qx^ 



{X = A,B), 



and the brackets without subscript (■ ■ ■) denote Gaussian averages, 

r°o dx 



(...) 



Analogously, the free energy per spin Q becomes 



f = — (1 - U) (1 - Qb) — 



(18) 



(19) 
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- ^ (In 2 cosh Ha) - ^ (In 2 cosh Hb) . 



1 - U? + (1 - gs?] + ^ruAms + | + ml) 



(20) 



III. THE STABILITY OF REPLICA-SYMMETRIC SOLUTION 



The validity of the RS solution (fTTj) rests on the applicability of Laplace method used to 
perform the integrations over m and q variables for large A^. The integral converges only 
if the stationary point (jl3j) is a minimum of 0, i.e., only if the eigenvalues of the Hessian 
matrix formed by the second derivatives of the function given by equation (fTUj) with respect 
to the m and q variables are all positive. We can equivalently consider as a function of 
A variables, related to m and q variables by means of Eq. ()12j). We will follow the latter 
approach because it leads to simpler calculations. The Hessian is a [n(n+ 1)/2] x [n(n + l)/2] 
matrix whose elements are 2x2 matrices given by 



a/3 



G 



where 



^AA ^AB 
\ ^BA ^BB 



Q«(/37) 



^AA ^AB 



\ ^BA 



^ Q(a/3)(7'5) 



BB 



^AA ^AB 

^(a/3)(7<5) ^(q/3){75) 
\ ^BA ^BB 



(21) 



G 



XY 



dX^dXy 



(/37)' 



G 



(a/3)(75) 
XY 



iX,Y = A,B). (22) 



At the stationary point of the RS solution |T7jl there are seven different types of 2 x 2 
elements of the Hessian matrix. We denote these elements by 



= A, 

Q(a/3)7 ^ g 



= B, 



Qa{al3) ^ Q 



Q{a/3)a ^ Q 



Qa{/37) ^ 13 



Q{a/3)(a7) ^ Q{al3)(yS) = 



(23) 



where the indices a, /3, 7 and 6 are all distinct and the tilde denotes the transpose of the 
matrix. We do not quote the lengthy expressions for these elements because only their linear 
combinations are needed in the calculation of the eigenvalues. 

The eigenvalues of the Hessian matrix can now be determined by finding the eigenvectors 
that divide the space into orthogonal subspaces closed to the permutation operation. The 
procedure are analogous to the case of the SK model j5| except that now the elements of 
the Hessian matrix are 2x2 matrices ((221) • These eigenvectors are 2^: n{n — 3) transversal 
or replicon eigenvectors depending on two replica indices, 4(n — 1) anomalous eigenvalues 
depending on a single replica index, and 4 longitudinal eigenvectors independent of replica 
indices. 
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The eigenvalues associated with the transversal eigenvectors are found to be the eigen- 
values of the 2x2 matrix 

T = P - 2Q + R, (24) 

with elements 

Tn = (1 - 2qA + ta) - {l3J'f{l - 2qA + r^)', (25) 
Ti2 = = -if3J)\l - 2qA + r^)(l - 2qB + r^), (26) 
T22 = (1 - 2qB + rs) - iPJ'fil - 2qB + rsf, (27) 

where 

tx = (temh' Hx) , rx = (temh' Hx) , {X = A,B). (28) 
The necessary and sufficient condition for all the eigenvalues to be positive are 

Tn + > and TnTsa - T^2 > 0, (29) 
which are equivalent to the conditions 

Ti = 2 - if3J')\l - 2qA + ta) - (/9J')'(1 - 2gB + Tb) > 0, (30) 
T2 = [l-{(3J')\l-2qA + rAm-{f3J'y{l-2qB + rB)] 

- i(3J)\l - 2qA + rA)(l - 2qB + Tb) > 0, (31) 



in agreement with previous studies 10lll3|. A RS solution satisfying these conditions will 
be called transversally (T) stable, and T unstable otherwise. 

The eigenvalues associated with anomalous and longitudinal eigenvectors are the same in 



the limit n 0. They are found to be the eigenvalues of the 4x4 matrix 

A-B D-C 

2C - 2D P - 4Q + 3R 

where 



(32) 



Lu = (1 - Qa) - PJ'^il - qA? + 2(/3J')'(tA - mA)\ (33) 

1^22 = (1 - qs) - I3J'q{1 - qs? + 2{l3J'f{tB - mB)\ (34) 

L12 = L21 = -/?Jo(l - - qs) + 2(/3J)2(t^ - mA){tB - m^), (35) 

1^13 = = {tA - mA)[l - PJoil - qA) - - Ma + Br^)], (36) 



1 



^L42 = {tB - mB)[l - PJoil - qs) - {^J')\\ - Mb + Sr^)], 



>32 



-PMtA - m^)(l - qb) - {l3Jf{tB - mB){l - Ma + Sr^), 



L2A 

-^23 
L44 

i^34 = i^43 = - Ma + 3r^)(l - 4gs + Sr^) + 2/3Jo(tA - mA){tB - ms). (42) 

The characteristic equation has the form 



;i - AqA + 3rA)[l - {(3J')\l - Ma + Sr^)] + 2/3J^(tA - niAf, 
;i - 4gB + 3rB)[l - - 4gB + Sr^)] + 2/3J^(tB - m^)', 



(37) 

(38) 

(39) 
(40) 
(41) 



A"^ — aiA'^ + a2A^ — asA + 04 = 0, 



(43) 



where the coefficients a„ are n-th order traces of the matrix L. A numerical study of 
equation ()43p shows that the eigenvalues are complex for some values of the parameters of 
the model, in contrast with one-sublattice SK model in which the anomalous and longitudinal 
eigenvalues neve, beccne eo.plex g Even thongh the Hessian n,atnx CD « > 1 .s 
real and symmetric, in the limit n — > there is no guarantee that the eigenvalues will 
be real. In fact, complex longitudinal and anomalous eigenvalues also arise in the spin 
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In 



1 one-sublattice infinite-range spin-glass model with crystal-field anisotropy 
general, therefore, the stability condition should require the real part of the eigenvalues to 



be positive. According to the Hurwitz criterion 



the necessary and sufficient condition 



for all the roots of equation 
Di = ai> 0, D2 



to have positive real parts are 
aia2 — 03 > 0, 

1 a2 



(44) 



Oi 03 
1 02 04 
ai 03 



03/^2 — ala^ > 0, 



Da 





03 








1 


02 


04 








ai 


03 








1 


02 


04 



a^Ds > 0. (45) 



These condition are equivalent to the following four conditions: 

Li = oi > 0, 
L2 = D2 = aia2 - 03 > 0, 

= Ds 
L4 = 04 > 0. 



010203 — O3 — 040^ > 0, 



(46) 
(47) 
(48) 
(49) 
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A RS solution satisfying these conditions will be called longitudinally (L) stable, and L 
unstable otherwise. 

IV. RESULTS OF THE STABILITY ANALYSIS 

In this section we present the results of the stability analysis of the RS solution for 
different values of the parameters of the model. Since the Hamiltonian is invariant under 
the simultaneous transformations 

H ^ -H, Si ^ -Si, (50) 

it is sufficient to consider fields H > 0. For H = only one of the two solutions related by 
the global inversion symmetry has to be considered. 

A. Zero applied field 

1. Ferromagnetic intersublattice interaction 

In zero applied field {H = 0) and ferromagnetic intersublattice interactions (Jq > 0) the 
solutions of the set of equations (fTTj) are of the form 

mA = mB = m, = = q. (51) 

Three types of solutions are possible: 

• Paramagnetic (P) solution: m = 0, g = 0. 

• Spin Glass (SG) solution: g > 0,m = 0. 

• Ferromagnetic (F) solution: g > 0, m > 0. 

Fig. ^ shows the lines delimiting the regions where different types of solutions can be found 
in the plane of temperature versus Jq + Jg- 

The P solution is always possible. However it is L stable only above the line (b) and the 
left portion of line (a), and T stable above line (a). The L instability of P solution occurs 
due to the violation of the condition (021), which is given in the case of P solution by 

L4 = [1 - P\J" - J'm - P'{J" + J')][l - P{4 - JoW - P{4 + Jo)] > 0. (52) 
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For (Jo + J'q)/\/ J'^ + J'^ < 1/2, the second factor in becomes negative below line (a). 
Thus the left portion of line (a) is determined by 

/32(J'2 + J2) = 1. (53) 

On the other hand, for (Jq + Jq) / \^J^~+~J^ > 1/2 the fourth factor in becomes negative 
below line (b). Thus the equation for line (b) is 

/?(J^ + Jo) = l. (54) 

The T instability of the P solution is due to the violation of the condition (jHlj) . which is 
given in the case of P solution by 

The second factor in becomes negative below line (a) for all values of Jq + J'q. Thus 
line (a) is given by equation (jK!^ for all Jo + Jq- The T and L instabilities of the P solution 
occur simultaneously on the line (a) for (Jq + Jq) /\/ J"^ + J'^ < 1/2. 

The SG solution is possible only below line (a). It is T unstable throughout this region 
and L stable to the left of line (c). The L instability of the SG solution occurs due to the 
violation of the condition (j^9|) . which is given in the case of SG solution by 

L4 = (1 - - 4g + 3r)2 [1 - - J^){1 - 4q + 3r)][l - + J^){1 - 4q + 3r)] 

X [1 - /3(J^ - Jo)(l - g)][l - P{Jo + Jo){l - q)] > 0. (56) 

For (Jo + Jo)/-\/J^ + J'^ > 1/2 the last factor in (fS^ becomes negative to the left of line 
(c). Thus the equation determining line (c) is 

/5(J^ + Jo)(l-g) = l. (57) 

The F solution is possible only between lines (b) and (c). It is L stable throughout this 
region but T stable only above line (d). The T instability of the F solution occurs due to 
the violation of the condition (jHT|) . which is given in the case of F solution by 

T2= [1-/52(J'2- J2)(i_2g + 3r)][l-/?2(j'2 + j2)(i_2g + 3^)] >0. (58) 

For (Jo + Jq)/ VJ^+J^ > 1/2 the second factor in becomes negative below line (d). 
Thus line (d) is described by equation 

/32(J'2 + J2)(l-2g + 3r) = 1. (59) 
10 



Rejecting solutions that are L unstable, we conclude that the P phase is located above 
lines (b) and left portion of line (a), the SG phase between the left portion of line (a) 
and line (c), and finally the F phase between lines (b) and (c). The SG phase, and the F 
phase between lines (c) and (d), are T unstable, indicating the need for a replica-symmetry- 
breaking solution in this region. The transition line (c) will change to a vertical line if such 



a solution is considered 



24| . We mention that, as should be expected, in the case Jo = 



and J = 0, or J/i = and J' = 0, these results reduce to those of one-sublattice SK model 



2. Antiferromagnetic intersublattice interaction 

The Hamiltonian in zero applied field {H = 0) is invariant under simultaneous trans- 
formations 

Jo-^-Jo, Si-^-Si (lEB). (60) 

In fact, we can check explicitly that all the expressions for the RS solution, including those 
of stability conditions, are invariant under simultaneous transformations 

Jo — ^ -Jo, mB — > -TUB- (61) 

Thus the case of antiferromagnetic intersublattice interaction Jq < is completely equivalent 
to the case of ferromagnetic intersublattice interaction — Jq > by replacing by —mB- 
This means that the F solution is replaced by the antiferromagnetic (AF) solution 

niA = -mB = m, Qa = Qb = q- (62) 

The results displayed in Fig. [T] remains valid, with Jq replaced by — Jq and F solution by 
AF solution. 



B. Non-zero applied field 

1. Ferromagnetic intersublattice interaction 

In non-zero applied field {H > 0) and ferromagnetic intersublattice interactions (Jq > 0), 
only the paramagnetic (P) solution is possible for the set of equations ()17p. which are of the 
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form 

rriA = trb = m > 0, = ?b = ? > 0. (63) 

This solution is always L stable, but becomes T unstable for low temperatures due to the 
violation of the condition (|31|). which in this case is also given by Eq. (j58p . The instability 
line is given by Eq. (jSHI), illustrated in Fig. |21 for the case J' / J = 1, Jq/ Jq = 1/2 and 
(Jo + Jo)/V + J'^ = 2 . As should be expected, in the case Jq = and J = 0, or 
Jq = and J' = 0, these results reduce to the de Almeida-Thouless line of one-sublattice 
SK modelQ. 

2. Antiferromagnetic intersublattice interaction 

In non-zero applied field {H > 0) and antiferromagnetic intersublattice interactions {Jq < 
0), two types of solutions to the set of equations (fTTj) are possible: 

• Paramagnetic (P) solution: rriA = rriB = m > 0,qA = Qb = Q > 0- 

• Antiferromagnetic (AF) solution: 7^ fnB,qA 7^ Qb- 

For (—Jo + Jq)/V .P + J'"^ < 1/2 only P solution is possible. This solution is always 
L stable but becomes T unstable at low temperatures due to the to the violation of the 
condition (jHTj) . which in this case it is also given by ()58|) . The instability line is given by 
Eq. dMD. 

For (— Jq + Jq)/V J2 + J'2 > 1/2, AF solution also becomes possible. Fig. El shows the 
lines delimiting the regions of existence and stability of each type of solution for — Jq/Jq = 
1/2. The P solution becomes L unstable below line (a) due to the violation of the condition 
In the case of P solution this condition is given by 

L4 = [(l-g)(l-4g + 3r)+2(t-m)2]2{[l-/5(J^ + Jo)(l-g)][l-/?(j'2_ J2) 
x(l - 4g + 3r)] + 20%^ + Jo){J" - J'){t - m)'} {[1 - /5(J^ - Jo)(l - q)] 
X [1 - /5( J'2 + J2)(l - 4g + 3r)] + 2/33( J^ - J,)(J'^ + J^)(t - mf] > 0. (64) 

The first factor in (jMj) becomes negative inside line (a). Therefore the equation determining 
line (a) is 

[1-/5(J^- Jo)(l-g)][l-/3(J'^- J^)(l-4g + 3r)] + 2/5Vo-'/o)(J"- J')(t-m)2 = 0, (65) 

12 



which is in agreement with previous study '125. The P solution is T unstable below line (b). 
This instability occurs due to the violation of condition (jHTj) . given in this case by Eq. (j^ . 
caused by the second factor. Therefore the line (b) is determined by equation (|59|). The AF 
solution is possible only inside line (a). It is L stable throughout this region and T unstable 
below line (c). This instability is due to the violation of condition (j3ip . Therefore line (c) 
obeys the equation 

[l-{(3J')\l-2qA + rA)][l-{f3J'f{l-2qB + rB)] = {(3J)\l-2qA + rA){l-2qB + rB). (66) 

Rejecting solutions that are L unstable, we conclude that P phase exists outside and AF 
phase inside line (a). The P solution becomes T unstable below line (b) and AF solution 
below line (c), which meet smoothly on the line (a). In the region below lines (b) and (c) it 
is necessary to consider replica symmetry breaking solution, which will presumably change 
line (a) in this region. 

For sufficiently large values of the ratio —Jq/Jq the model can exhibit first-order transition 
from AF phase to the P phase 12]. As an example, we consider the case —Jq/Jq = 5 shown 
in Fig. m The P solution is L stable inside line (a) given by Eq. ()65|). and T stable above line 

(b) given by Eq. (jUBj) . There is one AF solution inside line (a) and two distinct AF solutions 
between lines (a) and (d), as illustrated in Figs. El^a) andEfb) for ksT/ ^/ J"^ + J'^ = 1. One 
of the AF solutions, corresponding to dotted lines in Fig. is L unstable due to the violation 
of the condition (jl^ . as shown in Fig. Efc). The transition between AF phase and P phase 
is first order, determined by equating the free energies of L stable AF and P phases, as 
shown in Fig. Efd). The first-order transition line is shown as dotted line in Fig. which 
ends at the tricritical point TCP. The L stable AF solution becomes T unstable below line 

(c) due to the violation of condition (jHH), and it is determined by Eq. We conclude 
that in Fig. |3]the P phase exists outside and AF phase inside lines (a) and (e), which meet 
smoothly at the tricritical point TCP. The P phase becomes T unstable below line (b), and 
AF phase below line (c). Notice that the lines (b) and (c) are discontinuous across first-order 
transition It is likely that the first-order transition line will change in this part of phase 
diagram once the replica-symmetry-breaking solutions are considered for P and AF phases. 
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V. CONCLUSIONS 



In this paper we have investigated the stabihty of the RS symmetric solution of the 
two-sublattice generahzation of the SK infinite-range spin-glass model. We have derived 
stability conditions for transversal fluctuations in agreement with previous investigations, 
and we have extended previous study of the stability against longitudinal or anomalous 
fluctuations. The eigenvalues associated with such perturbations are in general complex. 
We generalized the usual stability condition by requiring the real part of these eigenvalues 
to be positive. The necessary and sufficient stability conditions were found using the Hurwitz 
criterion for all the roots of the secular equation to have positive real parts. These conditions 
allowed us to select one RS solution among those that are transversally stable. We believe 
that the generahzed stabihty condition should also be useful in other spin glass models where 
eigenvalues associated with longitudinal and anomalous perturbations become complex. 
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FIG. 1: Regions of the zero-field phase diagram where different types of solution are possible. For 
ferromagnetic intersublattice interaction (Jq > 0), the P solution is L stable above line (b) and the 
left portion of line (a), the F solution between lines (b) and (c), and SG solution between the left 
portion of line (a) and line (c). The P solution is T stable above line (a), the F solution between 
lines (b) and (d) . The SG solution is T unstable between the left side of line (a) and line (c) , and F 
solution between lines (c) and (d). For antiferromagnetic intersublattice interaction (Jq < 0) the F 
solution is replaced by AF solution and the label of horizontal axis by — Jq + Jq. The temperature 
and energy units in the axis are such that VJ^ + J'^ = 1 and A;^ = 1 
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H 

FIG. 2: Region of stability in the presence of a field for the case of ferromagnetic intersublattice 
interaction (Jq < 0). The values of parameters are J' / J = 1, Jq/Jq = 1 and ( Jq + Jq) I\/ J"^ + J'"^ = 
1/2. There is only the P solution which is always L stable but becomes T unstable below the de 
Almeida-Thouless line AT. The temperature and field units in the axis are such that \/J2TJ^ = 1 
and A;b = 1 
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FIG. 3: Regions of stability and existence of different solutions in the temperature versus field 
phase diagram for the case of antiferromagnetic intersublattice interaction (Jq < 0). The values 
of parameters are J' / J = 1, -J'q/Jq = 1/2 and (- Jq + Jq)H P + J'^ = 2. The P solution is L 
stable outside line (a) and T stable above line (b). The AF solution is possible only inside line (a) 
and it is always L stable, but becomes T unstable below line (c). The temperature and field units 
in the axis are such that \/ J"^ + J'"^ = 1 and ks = l 
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FIG. 4: Regions of stability and existence of different solutions in the temperature versus field 
phase diagram for the case of antiferromagnetic intersublattice interaction (Jq < 0). The values of 
parameters are J' / J = 1, —J^/Jq = 5 and (—Jo + Jq)/V J"^ + J'^ = 2. The P solution is L stable 
outside the line (a) and T stable above line (b). There is one AF solution inside line (a) which 
is always L stable, and two distinct AF solutions between lines (a) and (d), one L stable and the 
other L unstable. The L stable AF solution is also T stable above line (c). The line (e) is the 
first-order transition line and TCP is the tricritical point. The temperature and field units in the 
axis are such that V J'^ + J'^ = 1 and A;^ = 1. 
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FIG. 5: Field behavior of various quantities for the fixed value of temperature k^T j^J + J'^ = 1 
in the phase diagram of Fig. ^ In (a) and (b) the order parameters of the P and AF solution 
are shown, in (c) the L stability condition L4 of the AF solution, and finally in (d) the free-energy 
per spin of P and AF solutions. The L unstable AF solution are represented by dotted lines, and 
corresponds to the upper portion of the van der Waals loop in (d). The first-order transition is 
determined by the intersection of L stable AF and P solutions, as depicted in (d). The temperature 
and energy units in the axis are such that \/ + J' ^ = 1 and ks = 1- 
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